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ABSTRACT 

We present the results of a campaign to monitor active M dwarfs using infrared spectroscopy, 
supplemented with optical photometry and spectroscopy. We detected 16 flares during nearly 50 
hours of observations on EV Lac, AD Leo, YZ CMi, and VB8. The three most energetic flares also 
showed infrared emission, including the first reported detections of P/3, P7, He I A10830A and B1-7 
during an M dwarf flare. The strongest flare (Au = 4.02 on EV Lac) showed emission from H7, 
H<5, He I A4471A, and Ca II K in the UV/blue and P/3, P7, PS, Br 7 , and He I A10830A in the 
infrared. The weaker flares (Am = 1.68 on EV Lac and AU = 1.38 on YZ CMi) were only observed 
with photometry and infrared spectroscopy; both showed emission from P/3, P7, and He I A10830A. 
The strongest infrared emission line, P/3, occurred in the active mid-M dwarfs with a duty cycle of 
~3-4%. To examine the most energetic flare, we used the static NLTE radiative transfer code RH 
to produce model spectra based on a suite of one-dimensional model atmospheres. Using a hotter 
chromosphere than previous one-dimensional atmospheric models, we obtain line ratios that match 
most of the observed emission lines. 

Subject headings: line: formation - stars: activity - stars: chromospheres - stars: flare - stars: late- 
type 



1. INTRODUCTION 

M dwarfs are notorious for dramatic flares, presumably 
caused by magnetic reconnection in their atmospheres. 
Weaker analogs of these flares are present on the Sun, 
where the surface magnetic field is thought to be pow- 
ered by the rotationally induced shear b etween the ra - 
diative and convective layers of the Sun ([Parked fl955l ). 
For f ully convective M dwarfs (spectral types M3 and 
later: IChabrier fc Baraffel fT997'). strong magnetic fields 
are created and sustained sole ly through turbu lence and 
rotation within the star (e.g.. iBrowningl 12008). The dy- 
namo powering the magnetic fields which result in flares 
on these M dwarfs is not fully understood, but additional 
observations constraining the chromospheric heating are 
essential to a coherent picture relating stellar magnetic 
fields to flares. 

The atmospheric heating during a flare results in emis- 
sion from many wavelength regim es, and flares have been 
well-observed in the X-ray (e.g.. lOsten et al J 1201 Oft, u l- 
traviolet (e. g..lRobinson et al.ll2005t lHawlev et 81.11 2007), 
optical le.g. jKowalski et al.l2009tlWalkowicz et al.ll2011ft 
and radio (e.g.. iStepanov et al.l 120011 : lOsten fc Bastianl 
2008) . The combination of observations at these different 
wavelengths, especially when obtained as part of multi- 
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wavelength flare monitoring campaigns, has informed our 
interpretation of the p hysics underlying these dramatic 
emiss ion events (e.g., lHawlev et al.l 120031 : lOsten et al.l 
12005ft . However, to date there have been no concerted 
efforts to observe infrared emission lines from flaring 
stars. Quiescent M dwarfs are particularly bright in the 
near-infrared portion of the spectrum, so emission from 
the hydrogen Paschen and Brackett series and the He 
I A10830A transition are both easily observable and es- 
sential to probing different atmospheric heights. These 
infrared emission lines lines are particul arly useful for ex- 
amining accretion in T Tauri stars (e.g. lBarv et alJl2008t 
IVaccafe SandlllMl . 

In quiescent (not flaring) active M stars, high reso- 
lution (R > 20,000) spectra have shown weak ab sorp- 
tion from P/3 in AU Mic ([Short fc Dovld Il998bft : Pe 
was als o seen in absorption in six out of ten active M 
dwarfs ([Houdebine et all 12009ft . Emission from higher- 
order Paschen lines has only been detected in a few 
serendipitous observations at the far red end of op- 
tical spectra. The first occurred during a survey to 
classify photometric ally selected late-M and L dwarfs. 
iLiebert et al.l (|1999ft observed Paschen emission (PS - 
Pll) between 8800A and 10500A in a R-4300 spectrum 
of the M9.5 dwarf 2MASSW J0149090+295613. The 
flare also showed a variety of optical emi ssion lines, but 
had n o evidence of continuum emission. I Schmidt et al.l 
(2007) subsequently observed Paschen emission lines in 
a R~2000 spectrum of the M7 dwarf 2MASS J1028404- 
143843. The flaring spectrum included strong contin- 
uum enhancement of the entire spectrum blueward of 
9200A in addition to many emission lines. PS - Pll were 
again identified, with equivalent widths (EW) of 2-5A. 
For both of these observations, there was no correspond- 
ing photometry, which prohibited the characterization of 
the overall strength and duration of the flare. 
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iFuhrmeister et all (|2008l) observed PS - Pll on the 
M5.5 dwarf CN Leo during a large-amplitude flare with a 
total duration of about 45 minutes. Their data included 
R/^40,000 spectra over the range 3000A to 10500A. Line 
strengths were not given for the Paschen lines, but 
inspection of the 3 consecutive 100 s exposures near 
the peak of the f l are s hows a decay in line strength. 
IFuhrmeister et all ([2010D used one-dimensional atmo- 
sphere models to examine the emission from the flare, 
finding that a single model can reproduce most, but not 
all, flare emission lines. 

An unexplored region of the spectrum during M dwarf 
flares, both in observations and modeling, is the 1.0 to 
2.5/im range, which contains the lower order Paschen 
lines, higher-order Brackett and Pfund lines, and He I 
A10830A. We report on the results of our campaign to 
observe active M dwarfs in this wavelength regime with 
simultaneous photometric monitoring. The data include 
3 flares with infrared line emission; observations of the 
strongest flare also include blue optical spectra. Using 
these data, we quantify the duty cycle for infrared flare 
line emission and examine the relative line strengths dur- 
ing the evolution of each of th e three flares. Congru- 
ent with previou s stu dies (e.g.. lHawlev fc Fisherl 119921; 
iWalkowicz et all 120081 : IFuhrmeister et al.ll2010( )~we use 
one-dimensional atmosp heres and the sta tic NLTE ra- 
diative transfer code RH (Uitcnbrock 200 Q) to model the 
line flux ratios of the largest flare. 

In Section [5] we discuss our targets and observations 
and in Section [3] we describe our methods for flare iden- 
tification. Individual flares are examined in Section [4] 
together with our duty cycle estimate. Empirical atmo- 
spheric models that produce infrared line emission are 
discussed in Section [5] 

2. OBSERVATIONS 

Our targets include three well-known mid-M flare 
stars; AD Leo, EV Lac, and YZ CMi, in addition to one 
active late-M dwarf, VB 8. Magnitudes, coordinates, and 
the duration of our observations for each target are given 
in Table [1] To both detect IR emission and characterize 
a flare, we need at least infrared spectra and one band of 
optical photometry, but some nights include additional 
data. The targets, times observed, and instruments used 
each night are detailed in Table [5] 

2.1. Flare-cam on the ARCSAT 0.5-m 

Photometry for all of the nights was obtained us- 
ing Flare-cam on the Astrophysical Research Consor- 
tium Small Aperture Telescope (ARCSAT). ARCSAT 
was formerly used as the photometric calibrating tele- 
scope for the Sloan Digital Sky Survey ([York et al.ll2000l : 
I Tucker eT al. 2006). Its location at Apache Point Obser- 
vatory (APO) makes it an ideal telescope for obtaining 
simultaneous data with the ARC 3.5-m telescope, which 
we used for our infrared spectroscopy (see Section |2~3|) . 
Flare-cam is equipped with ugri filters, and the CCD is 
optimized for observing flare st ars because of its good 
blue response and fast readout (Hilton 201l|). Exposure 
times and filters are given in Table [5] 

The data were reduced using standard IRAF routines 
combined with a custom python code that tracked the 
change in each star's position over the course of the night. 



The magnitudes were calibrated using differential pho- 
tom etry wi t h resp ect to the brightest stars in the image. 
See IHiltonl (l20Tlh for more details on the photometric 
reductions. 

2.2. NMSU 1-m 

For four of the eleven nights of observations listed in 
Table [2] we also obtained £/-band photometry using the 
NMSU 1-m, a robotically op erated telescope located at 
APO ([Bolt zman et"aLl [2010). Typical exposure times 
were 4 -10 seconds, and readout was 10 seconds. The 
reductions were performed using an automated pipeline 
which measured the magnitude of the flare star with re- 
spect to several background stars. 

2.3. TripleSpec on the ARC 3.5-m 

Infrared spectra were obtained with the TripleSpec in- 
strument on the ARC 3.5-m telescope at APO. Triple- 
Spec is a cross-dispersed ne ar infrared spectrog raph that 
covers 0.95 to 2.45 microns ([Wilson et al.ll2004l) . We used 
the 1.1" slit, resulting in a resolution of R~3500. All 
data were obtained using an A/B nod pattern, shifting 
the star along the slit in order to perform sky subtrac- 
tion. We obtained data for an AO calibrator star every 
^40 minutes in order to correct for the changing telluric 
absorption over the course of the night. Our typical ex- 
posure times were 2-4 seconds for AD Leo, EV Lac, and 
YZ CMi, and 30 seconds for VB 8. 

The data were reduced using a version of SpexTool 
modified to work w ith ARC 3.5-m TripleSpec data 
([Cushing et al.l [20041 ). We constructed telluric correc- 
tion spectra from our AO standards usi ng the routine 
included in SpexTool ([Vacca et al.l 12003]), but modified 
the remaining post-processing routines to automatically 
process each spectrum instead of using the provided GUI 
interface. Although the formal residuals of our wave- 
length solution were 0.5-1 pixel, the curved, tilted orders 
of the spectra impose additional systematic effects in the 
wavelength calibration. 

We detected P/3, PS, P7, Br 7 , and He I A10830A dur- 
ing the most energetic flare observed. No higher order 
Brackett or Pfund emission was detected in any of our 
spectra. We measured the equivalent widths (EW) us- 
ing regions defined individually for each line in order 
to include all observed flux; for P/3, PS, P7, and He 
I A10830A these were 10-20A wide (6-12 pixels; 0.001- 
0.002 (im) and for P>rj the line region was 40A wide (14 
pixels; 0.004 //m). Continuum regions were defined as 
iO.Ol^m on either side of each line. Quiet and flare pro- 
files for the five lines are shown with the regions used for 
line measurements in Figure [TJ The spectrum surround- 
ing and underlying each of the emission lines contains 
many other molecular and atomic features so the EWs of 
the emission lines are not zero even in quiescence. The 
EW measured in the quiescent spectrum is subtracted 
from each flare measurement. 

We could not determine absolute line fluxes directly 
from individual TripleSpec spectra because the observed 
flux in each spectrum varies due to the movement of 
the star on and off the slit during the nod pattern. We 
used a method similar to the x factor of lWalkowicz et all 
( 2004) to calculate absolute line flux by multiplying mea- 
sured equivalent widths (which do not depend on the 
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TABLE 1 

Flare Star Observations 



Name 


ST 


J 


K 


tobs (h:m) 


Nflarcs 


tflarc (h:m) 


Nflarcs IR 


tflare IR (h:m) 


Frac m 


YZ CMi 


M4.5 


6.58±0.02 


5.70±0.02 


16:36 


8 


1:55 


1 


0:08 


0.008 


AD Leo 


M3 


5.45±0.02 


4.59±0.02 


12:18 


1 


0:31 





0:00 





EV Lac 


M3.5 


6.11±0.03 


5.30±0.02 


15:21 


6 


3:56 


2 


1:14 


0.081 


VB 8 


M7 


9.78±0.03 


8.82±0.02 


4:37 


1 


0:15 





0:00 





M3-M4.5 








44:15 


15 


6:22 


3 


1.37 


0.031 


total 








48:52 


16 


6:38 


3 


1.37 


0.028 



Note. — The total time observed is given in column 5 
(tfiare), and the time with observed infrared line emission 
gives the fraction of time each object spent with infrared 



(t bs), the total time each object spent in flare is given in column 7 
is given in column 9 (t b s ; see Figurcs[3] \E\ and[6j. The last column 
line emission. 



TABLE 2 
List of Observations 



UT date 


Target 


Telescope/Instr . 


Filter 1 


ET (s) 


UT Time (h:m) 


Time (h:m) 


tobs (h:m) 


Nflarea 


2009 Feb 4 


YZ CMi 


APO TripleSpec 2 






4 


3:20 to 7:10 


3:50 


1:42 


1 


2009 Feb 4 


YZ CMi 


ARCSAT flare-cam 


u, g 




7, 2 


5:30 to 7:12 


1:42 






2009 Feb 4 


AD Leo 


APO TripleSpec 






2 


7:14 to 13:46 


6:32 


5:20 


1 


2009 Feb 4 


AD Leo 


ARCSAT flare-cam 


g 




5, 1 


7:18 to 12:38 


5:20 






2009 Apr 14 


AD Leo 


APO TripleSpec 






2-5 


2:54 to 7:44 


4:50 


1:18 





2009 Apr 14 


AD Leo 


ARCSAT flare-cam 


g 




10-100 


2:28 to 4:12 


1:44 






2009 Oct 24 


EV Lac 


APO TripleSpec 






8 


1:41 to 6:51 


5:10 


3:35 


1 


2009 Oct 24 


EV Lac 


ARCSAT flare-cam 


u, g 




10, 1 


1:06 to 5:05 


3:59 






2009 Oct 24 


EV Lac 


DAO Spectrograph 3 






60-300 


2:26 to 9:05 


6:39 






2009 Oct 24 


EV Lac 


NMSU 1-m camera 


U 




4 


1:26 to 9:30 


8:04 






2009 Oct 27 


EV Lac 


APO TripleSpec 






8 


1:04 to 6:53 


5:49 


5:19 


3 


2009 Oct 27 


EV Lac 


ARCSAT flare-cam 


u, g 




5, 1 


1:34 to 9:13 


7:49 






2009 Oct 27 


EV Lac 


DAO Spectrograph 






200-420 


2:02 to 6:18 


4:16 






2009 Oct 27 


EV Lac 


NMSU 1-m camera 


U 




4 


1:12 to 9:18 


8:06 






2010 Apr 25 


VB 8 


APO TripleSpec 






60 


7:24 to 11:07 


3:43 


3:12 





2010 Apr 25 


VB 8 


ARCSAT flare-cam 


g 




120 


7:57 to 11:06 


2:33 






2010 May 26 


VB 8 


APO TripleSpec 






30 


7:28 to 8:53 


1:25 


1:25 


1 


2010 May 26 


VB 8 


ARCSAT flare-cam 


g 




60 


6:22 to 8:54 


2:32 






2010 May 26 


EV Lac 


APO TripleSpec 






8 


9:01 to 11:39 


2:38 


1:50 


1 


2010 May 26 


EV Lac 


ARCSAT flare-cam 


u, g 




10, 2 


9:10 to 11:00 


1:50 






2010 Nov 27 


EV Lac 


APO TripleSpec 






5 


0:24 to 5:14 


4:50 


4:37 


1 


2010 Nov 27 


EV Lac 


ARCSAT flare-cam 


u, g. 


r 


5, 1, 1 


0:37 to 5:17 


4:40 






2011 Feb 14 


YZ CMi 


APO TripleSpec 






4 


2:07 to 6:08 


4:01 


4:01 


4 


2011 Feb 14 


YZ CMi 


ARCSAT flare-cam 


g> r, 


i 


1, 1, 1 


1:24 to 4:30 


3:06 






2011 Feb 14 


YZ CMi 


NMSU 1-m camera 


U 




10 


2:07 to 7:20 


5:13 






2011 Feb 15 


AD Leo 


APO TripleSpec 






2 


9:35 to 12:28 


2:53 


2:36 





2011 Feb 15 


AD Leo 


ARCSAT flare-cam 


g, r, 


i 


1, 1, 1 


9:26 to 10:36 


1:10 






2011 Feb 15 


AD Leo 


NMSU 1-m camera 


U 




4 


9:52 to 13:09 


3:17 






2011 Feb 15 


YZ CMi 


APO TripleSpec 






4 


1:16 to 3:08 


1:52 


1:07 


2 


2011 Feb 15 


YZ CMi 


ARCSAT flare-cam 


g> r, 


i 


1, 1, 1 


1:26 to 3:56 


2:30 






2011 Feb 15 


YZ CMi 


APO TripleSpec 






4 


4:49 to 9:34 


4:45 


4:45 




2011 Feb 15 


YZ CMi 


ARCSAT flare-cam 


g, r, 


i 


1, 1, 1 


4:50 to 5:27 


1:37 






2011 Feb 15 


YZ CMi 


NMSU 1-m camera 


U 




10 


1:59 to 9:41 


7:42 






2011 Feb 21 


AD Leo 


APO TripleSpec 






2 


8:41 to 12:16 


3:35 


3:04 





2011 Feb 21 


AD Leo 


ARCSAT flare-cam 


u, g. 


r 


5, 1, 1 


6:07 to 11:45 


5:38 






2011 Feb 22 


YZ CMi 


APO TripleSpec 






4 


2:08 to 7:09 


5:01 


5:01 


1 


2011 Feb 22 


YZ CMi 


ARCSAT flare-cam 


"> g: 


r, i 


6, 1, 1, 1 


1:49 to 7:40 


5:51 







Note. — Column 7 (Time) refers to the total duration of observations for each instrument, while column 8 (T b s ) gives the duration 
of observations that overlap for the photometry and infrared spectroscopy. Column 9 (Nfi ares ) gives the number of flares found from 
the photometry using the process described in Scction[3l 

The filter in bold was used to identify flares for those observations. 

2 R~3500; A ~ 0.95 to 2.45/-tm; sec section [Q] 

3 R~750; A ~ 3540 to 4710A; see section \2A\ 
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Fig. 1. — Spectra of EV Lac during quiescence (black) co-added 
from ~2 hours of exposure during UT 2009 October 24, and from 
8 co-added spectra with exposure times of 8 seconds each during 
the flare peak on UT 2009 October 27 (red). The emission lines 
are labelled and their central wavelengths are shown (dashed green 
lines). The regions used to measure the EW of each line are also 
shown (solid blue lines). The shifts of some lines from the nominal 
central wavelengths are likely due to small (1-2 pixel) systematic 
errors in the wavelength solution (see text). 



continuum level) by a calibrated continuum flux. Con- 
tinuum fluxes were obtained from a quiet, co-added, high 
S/N spectrum of each star normaliz ed to 2MASS pho- 
tometry using 2MAS S filter curves (jCohen et al.l 120031 : 
iSkrutskie et al.l[2006f ). While this method would not be 
feasible in the UV and optical due to white light emission, 
continuum enhance ment during flares is neg ligible in the 
JHKs passbands. IDavenport et al.l (|2011l) use a flare 
continuum model on an M3 star to predict that a flare 
with Am = 4 mags would produc e a A J < 10 millimags 
peak, and Tof flemire et al.l (|2011[) report no broad-band 
(J, H, Ks) continuum enhancements above a level of 5-8 
millimags during flares having similar t otal energy as re- 
orted her e. The variations detecte d bv lTofflemire et all 
201 If ) and IDavenport e t aL (20lT|) are smaller than the 
formal uncertainty quoted with the 2MASS magnitudes 
(20-30 millimags), so we assume that the variation be- 
tween the published 2MASS magnitudes and the magni- 
tudes of the M dwarfs during our observing was negligi- 
ble. 

2.4. DAO 1.8-m 

For two nights, we used the DAO 1.8-m telescope with 
the SITe5 CCD and spectrograph to observe EV Lac dur- 
ing a coordinated campaign with the telescopes at APO. 
Our setup resulted in a spectral resolution of R~750 and 
wavelength coverage from 3550A-4700A. We measured 
Ca II K, He I A447lA, and the Hydrogen Balmer series 
Hy and H<5. Exposure times for EV Lac ranged from 60 
to 420 seconds. Due to these relatively long integration 
times, additional cosmic ray cleaning was pe rformed with 
the LACOSMIC utility (jvan Dokkum|[200ll) . 

The spectra were wavelength-calibrated with a FeAr 
lamp and flux-calibrated using data from the standard 
star G191B2B, then spectrophotometrically calibrated 
by normalizing to the simultaneous U -band data. Equiv- 
alent widths are not useful for blue flare spectra because 
of the changes in the surrounding continuum flux during 
the flare. Instead, we measured absolute line fluxes di- 
rectly from the data. The values we use during the flare 
have the quiet line flux subtracted. 

3. IDENTIFYING FLARES 

Flares are most easily seen at blue and ultraviolet 
wavelengths, where the hot, white-light continuum emis- 
sion from the flare is in high contrast to the small 
amount of flux emitted from cool M dwarf ph otospheres 
(jLacv et al.l 119761: lHawlev fc Pettersenl 119911) . To iden- 
tify as many flares as possible, we used the bluest band 
of photometry available. This was typically u, but for 
some nights only U was available, and VB 8 was too 
faint to observe in U or u, so we used g-band data. The 
band used to identify flares for each set of observations 
is given in Table [5J 

Photometrically, flares are observed as excursions 
above the mean quiescent value of the star's flux, which 
can be any size or shape. Realistically, flare detection 
must take into account small variations in the contin- 
uum caused by observational effects and so a minimum 
duration and energy above the observed quiescent value 
is required. To identify indivi dual fl a res, w e used the 
custom IDL code discussed in Hilton] ([201 If) , which se- 
lects peaks that have at least three consecutive epochs 
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Fig. 2. — Peak magnitude in U-, u-, or g-band as a function of 
total flare energy in those bands. The flares from AD Leo (trian- 
gles), YZ CMi (squares), EV Lac (circles) and VB 8 (diamonds) 
observed in U (red), u (blue), and g (green) band are shown. The 
flares where IR line emission was observed are distinguished (solid 
symbols) and labelled with their dates. The flares with accom- 
panying IR line emission have the largest total flare energy and 
relatively high peak magnitudes. 



more than 3 standard deviations above the local quies- 
cent light curve. At least one of those epochs must be 
5 sigma above quiescence. We reviewed each flare by 
eye to confirm that the deviations from the mean were 
not caused by bad photometry. Over the course of 48.9 
hours of observations on four different stars, we observed 
a total of f6 flares, which are listed per star in Table [T] 
and per night in Table [2j Figure [2] shows the energy and 
peak magnitude of each flare. 

To identify flares which had associated IR line emis- 
sion, we examined the measured EWs of P/3 and He I 
AI0830A as a function of time during the flare. We found 
that the three most energetic flares, which occurred on 
EV Lac on UT 2009 October 27 and UT 2010 Novem- 
ber 27, and on YZ CMi on UT 2011 February 14, each 
showed infrared emission lines. These flares are discussed 
in detail in Section |4~T1 



4. CHARACTERIZING INFRARED FLARES 

4.1. Individual Flares 

2009 October 27 flare on EV Lac: We were observ- 
ing with all four instruments during the most energetic 
event, a Au = 4.023 magnitude flare on EV Lac on UT 
27 October 2009. The light curves for our observations 
are shown in Figure [3J The photometry (in U-, u-, and 
g-band) exhibits a typical flare light curve with a fast rise 
and exponential decay. The u-band flare emission lasted 
1.68 hours and released a total energy of 3.9 xlO 32 ergs. 

2 Although the Am represents a negative change in magnitude 
(corresponding with an increase in flux) we adopt a convention of 
Au = \u flare - u quiet \ throughout. 



The combination of optical (DAO) and infrared 
(TripleSpec) spectroscopy allows us to examine a total 
of nine emission lines - H 7 , US, He I A447lA, and Ca 
II K in the UV/blue part of the spectrum, and P/3, P7, 
PS, Br 7 , and He I A10830A in the infrared. Figure H 
shows the lightcurve of each emission line normalized to 
its value at t = 4.97 hours (the peak of u-band emission), 
and the ratio of each line to H7 for comparison of their 
evolution during the flare. 

The light curves for H 7 , Ed, and He I A447lA have 
a fast-rise exponential-decay shape similar to the pho- 
tometry. P7 and PS show a similar fast rise, but their 
decay is slower than the Balmer series lines. The P/3 and 
Ca II K emission both peak after the other Paschen and 
Balmer series lines, and exhibit an even slower decay af- 
ter their late peaks. Br 7 is similar to P/3 and Ca II K in 
its late peak, but seems to decay faster than any other 
line. This may be an observational effect, as it is by far 
the weakest line detected. Without a stronger detection, 
we assume that its ratio to the Paschen lines is constant 
throughout the flare. The He I A10830A emission shows 
a shape distinct from the rest of the lines - it remains 
nearly at its peak flux for 0.8 hours, approximately half 
of the duration of the flare in it-band. 

The slow decay during the gr adual phase is a well- 
known property of Ca II K (e.g. iBopp fc MoffettJ 11 973t 
lHawlev fe PettersenlH99ll: iFuhrmeister et al J 12008 ). but 
in this flare, He I A10830A emission traces a region that 
remains heated for an even longer portion of the gradual 
phase than Ca II K. This could be due to the Neupert ef- 
fect, where the line responds to the total cumulative flare 
heating for which the time integral of the U-band (white 
light emission) is often used as a proxy (|Hawlev et al.1 
Il995t lOsten et al.ll2005l) . Section [5] describes our efforts 
to model the emission lines from this flare. 

2010 November 27 Flare on EV Lac: We observed 
another flare with infrared line emission on EV Lac on 
UT 2010 November 27. The flare peaked at Am = 1.68, 
and over the course of t = 1.30 hours it emitted 5.5 x 10 31 
ergs in the u-band. We observed with both ARCSAT 
and TripleSpec during the flare, and have photometry in 
g and r-band in addition to the u-band data. The pho- 
tometry and the line flux lightcurves for P/3, P 7 , and He I 
A10830A are shown in Figure [5j There was no discernible 
emission in PS and Br 7 . 

This peculiarly-shaped flare contains three separate 
peaks in the u-band photometry. After the first and 
third peak, the flux seems to decay exponentially, but 
after the middle peak there is a gentle rise in the u-band 
flux. TripleSpec was taking observations of a standard 
star during the first peak of the flare, so it is unknown if 
the emission lines showed the same fast-rise exponential- 
decay as the first photometric peak. The rise in P/3 and 
P 7 line emission before and after the standard star gap 
suggests that those lines showed some emission between 
the first and second peaks of the flare. An observed in- 
crease in infrared line emission occurred ~0.2 hours after 
the second peak in the u-band photometry, tracing a gen- 
tle rise and decay. 

The shape of this flare is very different than that of the 
UT 2009 October 27 flare on EV Lac, and the relative 
line strengths are also different. In the previous flare, 
P/3, P 7 and He I A10830A emitted nearly the same peak 
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Fig. 3. — Data from the UT 2009 October 27 flare observed on EV Lac. Panels (a), (b), and (c) show U, u, and g-band light curves 
during the flare. The panels (d) - (1) show the variations in the optical and infrared lines. The flaring line flux (with a quiescent value 
subtracted) is shown along the left y axis, and for the infrared data equivalent width is shown along the right y axis. Panel (k) shows the 
time used as the flare duration in Table [T] (red line). The time of the peak -it-band flux is shown in each panel (blue dashed line). The 
infrared emission line fluxes are averaged over 8 exposures of 8 seconds each. The total exposure time for each point is 64 seconds, but 
including instrumental overhead and time to execute the nod pattern, the average cadence is 2.5 minutes. Small gaps in the infrared line 
data are due to standard star observations. 
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Fig. 4. — Line fluxes for the nine optical and infrared emission 
lines detected during the UT 2009 October 27 flare observed on EV 
Lac. In panel (a), the line fluxes are normalized to their flux at the 
time of peak emission for H7 (t = 4.97 hours). The normalization 
emphasizes the different decay patterns. Panel (b) shows the ratios 
of the other eight emission lines to the H7 flux. The IR lines were 
binned to the same time resolution as the optical lines. The vertical 
green dashed line shows the location of the flare peak in H7. 

flux. In this flare, He I A 10830 A peaked at twice the 
strength of the P/3 and P7 lines, indicating a different 
pattern of atmospheric heating during the two flares. 

2011 February 14 flare on YZ CMi: On UT 2011 
February 14, we observed a AU = 1.38 flare on YZ CMi 
with the NMSU 1-m, ARCSAT, and TripleSpec. The 
flare lasted for t = 0.5 hours and released a total U- 
band energy of 4.4 x 10 31 ergs. Figure [B] shows the U, 
g 7 and r-band light curves (i-band was also observed but 
showed no change during the flare) and line fluxes from 
P/3, P7 and He I A10803A. This is the lowest energy flare 
with any evidence of IR line emission, and the measured 
EW were small (0.05 to 0.2A), which provides a lower 
limit on the observability of IR line emission. With these 
small EW, it is difficult to compare the strengths of the 
emission lines; they are all the same strength within the 
uncertainties. 

The [/-band light curve shows a fast-rise exponential- 
decay shape with a precursor event 0.1 hours before the 
main peak. The IR emission does not show the precur- 
sor or the initial rise of the photometry. However, the 
co-added infrared measurements have an effective time- 
resolution of 2.5 minutes (due to the inclusion of time 
spent executing the nod pattern and readout), which is 
insufficient to resolve those features. 

4.2. How Often Does IR Line Emission Occur? 

In order to determine the expected rate, or duty cycle, 
of infrared line emission, we first defined detectable emis- 
sion as approximately la above the mean quiescent level. 
The length of time with detectable emission is shown for 
each flare as the red horizontal line on the P/3 light curve 
in Figures [3J [S] and [5] and given in Table [TJ The total 
time spent in emission for all three flares observed is 1.4 
hours (out of 48.9 possible hours), which corresponds to 
an IR flare emission duty cycle of 2.8%. Excluding VB 
8, we calculate an IR flare emission duty cycle of 3.1% 
(of 44.3 hours) for active mid-M dwarfs. 

We can also place a limit on the duty cycle using the 
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Fig. 5. — Data from the UT 2010 November 27 flare observed 
on EV Lac. Panels (a), (b), and (c) show the change in «-, g-, 
and r-bands during the flare, and panels (d), (e), and (f) show the 
flux (left axis) and EW (right axis) of the infrared emission lines 
detected during the flare. Panel (d) shows the time used as the 
flare duration in Table [T] (red line) . The three dashed blue lines 
show the times when the three flare peaks occurred in the ii-band. 
The infrared emission line fluxes are averaged over 8 exposures of 5 
seconds each. The total exposure time for each point is 40 seconds; 
the average cadence is 1.7 minutes. The gaps in the infrared line 
data are due to standard star observations. 
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in each panel (blue dashed line). The infrared emission line fluxes 
are averaged over 16 exposures of 4 seconds each. The total expo- 
sure time for each point is 64 seconds; the average cadence is 2.5 
minutes. The gap in the infrared line data is due to standard star 
observation. 



flare frequency distributions from iHiltonl f|2011h , which 
give the number of w-band flares per unit time for each 
flare energy. The u and U -band energies of the flares with 
accompanying infrared emission are all above 3 x 10 31 
ergs, corresponding to a flare frequency <0.1 per hour. 
Multiplying this emission time per flare by the flares per 
hour gives a duty cycle of <4.6%, in agreement with 
our independent estimate. A duty cycle of 2.8 - 4.6% 
represents an upper limit on detectable emission at this 
S/N and resolution, as our criterion requires only a small 
detection in the brightest line. 

5. ATMOSPHERIC STRUCTURE 

We used the st atic NLTE radiative transfer code RH 
(Uitenbroek 2001) to generate model spectra to compare 
with the emission lines observed in the UT 2009 October 
27 flare on EV Lac. We calculated model spectra based 
on one-dimensional atmospheres, using a 20-level hydro- 
gen atom, a 20-level calcium atom, and a 25-level he- 
lium atom. The multi-level atoms were required to gen- 
erate the lines observed, while the simplification to a one- 
dimensional atmosphere allowed us to examine a larger 
range of chromospheric structures without the compu- 
tationally intensive calculation s required by a d etailed 
treatment of flare physics (e.g.. lAllred et al.ll2006l) . 

For a starting atmosphere, we used a Nextgen photo- 
spheric model from a T = 3200K solar metallicity dwarf 
( Hauschil dt et al.lll999t ) and the corona o f the pr e-flare M 
d warf atmosphere model o f lAllred et al.1 (l2006f). S imilar 
to lHawlev fc Fisher! (IT992T ). iChristian et all (|2003h and 
iFuhrmeister et al.l (|2010D . we used chromospheres with 
a linear temperature rise in log column mass (log(col 
mass)) to connect the photosphere and corona. The lin- 
ear temperature rise is a simplification of the actual chro- 
mospheric structure during a flare, but it is useful for an 
initial investigation of the temperatures required to gen- 
erate emission lines at each atmospheric height. To pro- 
duce a suite of model atmospheres, we varied the column 
mass of the transition region (log(col mass)™), the col- 
umn mass of the temperature minimum region (log(col 
mass) Tmin) , and the temperature of the chromosphere 
at the bottom of the transition region (Ttr)- Figure [7] 
shows the temperature structure of a representative sub- 
set of the resulting atmospheres and illustrates the three 
quantities w e varied. 

Following IWalkowiczl (|2008f l. we adopted our initial 
ranges for Ttr, log(col mass)™, and log(col mass)Tmm 
from previous quiescent and flaring; M dwarf chromo- 
spher e models (IHawlev fe Fisher! 1 1992t iMauas fc Falchi 



1991 IHoudebine fe Stempeh lT997t IShort fe Dovle 



1998a! IWalkowicz et all 120081: Fuhrmci ster et al. 



Table [3] shows the range of parameters adopted for 
each of these three quantities, which differ from pre- 
vious parameter ranges only in Ttr- Previous model 
atmospheres have relatively constant Ttr ~ 10 000K, 
but our initial models with a range of Ttr = 10 000K 
to 20 000K underproduced Paschen emission relative 
to Balmer emission, and showed a trend of increasing 
Paschen emission with greater Ttr- We increased the 
temperature of our hottest models to Ttr = 30 000K 
in order to generate relatively more Paschen emission. 
Although Ttr extends to hotter temperatures, it is 
consistent with re sults from th e radia tive-hydrodynamic 
simulations from lAllred et alj ([2006), which show that 
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Fig. 7. — The atmospheric temperature structure as a function 
of log(col mass) for a subset of the atmospheres generated. Panels 
(a) and (b) show the entire atmosphere, while panels (c) and (d) 
focus on the chromosphere. The panels (a) and (c) show the models 
generated with a log(col mass)x m in — and panels (b) and (d) 
show log(col mass)y m i„ = — 1. Each panel shows 9 models, with 
a range of log(col mass)TR = —3.5 to —5.5 (blue dot-dashed lines) 
and a range of T TR = 20 000K to 30 000K (dotted green lines). 
Line fluxes ratios for the spectra generated from these models are 
shown in Figure [8] 



TABLE 3 

Model Atmosphere Parameters 



Parameter 


Range 


Best 


log(col mass)xij 


-5.5 to -3.5 


-5.5 to -4.5 


log(col mass)Tmin 


-3 to 





Tta 


10 000K to 30 000K 


25 000K to 30 000K 


T 

- 1 - m in 


2229K to 3264K 


3264K 



material at the base of the transition region can be 
heated up to T = 10 6 K. 

Comparing the strengths of the modeled lines to each 
other provides strong constraints on our suite of model 
atmospheres. The line flux ratios with respect to the H7 
line flux for the best models are shown compared to the 
median and range of observed line flux ratios in Figure 
In general, a deeper T m i„ (at log(col mass) = or — 1), 
a deeper transition region (at log(col mass) = —3.5 or 
—4.5), and a hotter chromosphere (with Ttr — 25 000K 
or 30 000K) better reproduce the line flux ratios observed 
during the flare. The line formation regions (where the 
contribution function for each line is greater than 25% of 
its peak value) for one model are shown in Figure [9] 

Most of the lines are produced over regions that in- 
clude the outer portion of the chromosphere, at log(col 
mass) = —3.5 and T = 25 000K. H7, the strongest emis- 
sion line we observed, is formed over the smallest portion 
of the chromosphere, with its highest temperature at T 
= 20 000K. He I A447lA and He I A10830A form over 
slightly different regions of the atmosphere, with He I 
A10830A tracing slightly higher temperatures. 

During the flare observations, the ratio of US to H7 
is relatively constant. Hi5 is slightly overproduced in the 
models compared to observations, but is similar in each 
of the models. He I A447lA emission is weak (with a 
ratio to H7 of ~ 0.1 to 0.3) in both the models and the 
observations. The ratios of the Paschen series lines and 
Brackett 7 to each other are relatively constant through- 



out the flare, and those ratios are well produced in every 
atmospheric structure. The ratio of the Paschen lines 
and Br7 to H7, however, is matched only in the mod- 
els with a log(col mass) Tmin = 0. Because these lines 
are sensitive to the hottest regions of the chromosphere, 
the advantage of the log(col mass)Tmin = is likely an 
increased amount of material at temperatures near T — 
20 000K due to a shallower slope in the chromosphere. 

The ratio of Ca II K to H7 and to the other hydro- 
gen series lines is best produced in the models with the 
deepest T m j„, a deep transition region (log(col masses 
= -3.5 or -4.5), and a hot T T r = 25 000K or 30 000K. In 
all other models, Ca II K is underproduced relative to 
the Paschen series lines. As shown in Figure [9l the Ca II 
K emission in the best-fit model is formed over a larger 
range of log(col mass) than any other line. The produc- 
tion of Ca II K in a region that includes the upper chro- 
mosphere is unusual compared to previous results; typi- 
cally, Ca II K emission during a flare is tho ught to last 
longe r because it is a lower temper ature line (Houdcbinc 
I2003t iCrespo-Chacon et all I2006D . Our cooler atmo- 
spheres, where Ca II K emission is formed only in lower 
temperature regions, do not produce enough Ca II K 
emission relative to Paschen series emission to match our 
observations. 

He I A10830A is underproduced in nearly every model. 
During the flare, its observed ratio compared to H7 in- 
creases from 0.1 to 0.5, while all our models show line 
flux ratios of 0.1 or less. This mismatch is apparently 
worse in the one of the other two flares observed; as dis- 
cussed in Section |H He I A 10830 A is stronger compared 
to P/3 and P7 (the two other lines observed) in the UT 
2010 November 27 flare on EV Lac. Simply raising the 
Ttr in our models produces too much Ca II K but no 
additional He I A10830A. In the Sun and similar stars, 
He I A10830A emission is produced in the upper chromo- 
sphere during flares as a result of he lium ionization via 
backwarming from coronal UV flux (Mauas et al.l [20051 : 
iSanz-Forcada fc Dupred l2008f ) . A similar process could 
be leading to the He I A 10830 A emission during M dwarf 
flares, but the details of backwarming from coronal emis- 
sion are not yet fully implemented in the RH atmosphere 
code. 

While one-dimensional atmosphere models can match 
the line flux ratios of most of the lines as a sequence of 
static snapshots, they cannot reproduce the time evo- 
lution of the flare. In our observations, Ca II K, the 
Paschen lines, Br7, and He I A 10380 A all rise relative 
to H7 during the decay phase of the flare. The best 
fitting models in our suite of atmospheres indicate an in- 
crease in Ca II K is always coupled with a decrease of 
the Paschen and Brackett lines. The time evolution of 
flares may involve different atmospheric components cov- 
ering the surface of the sta r with changing filling factors 
(e.g.. [Kowalski et al.ll2010T ). It is possible that a linear 
combination of two or three different one-dimensional at- 
mospheres with changing filling factors would reproduce 
the time-evolution of this flare. 

6. SUMMARY 

During nearly 50 hours of simultaneous photometric 
and spectroscopic observations on 4 active M dwarfs, we 
saw 16 total flares, 3 of them with accompanying infrared 
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Fig. 8. — Line flux ratio to H7 flux as a function of transition region temperature generated from the models shown in Figure [7] Each 
line flux ratio is shown in a different color and symbol, which is detailed at the top. The two rows show the line flux ratios from models 
with different log(col mass)T m i n . For comparison, the observed median and range of each line flux ratio is shown to the left of both rows. 
Each panel (a-f) shows models with a different log(col mass)-^ and three values of Ttr- The dotted lines connecting the model line fluxes 
are shown only to clarify the positions of crowded points. The models that best represent the data are found in panels (a) and (b), with 
T TR = 25 000K and 30 000K. 



emission lines. The strongest flare (Au = 4.02) occurred 
on EV Lac on UT 2009 October 27. It showed emis- 
sion from H 7 , Hi, He I A447lA, Ca II K, P/3, P 7 , PS, 
Br 7 , and He I A10830A. A weaker flare (Au = 1.68) on 
EV Lac on UT 2010 November 27 showed only emission 
from P/3, P 7 , and He I A10830A. Remarkably, the He I 
A10830A emission was twice as strong compared to P/3 
and P 7 as it was in the Am = 4.02 flare. The weakest flare 
with infrared emission (AU = 1.38) occurred on YZ CMi 
on UT 2011 February 14; P/3, P 7 , and He I A10830A were 
just above their detection limits. We estimate a duty cy- 
cle of 2.8% to 4.6% for observing the strongest infrared 
emission line (P/3) during flares on active mid-M dwarfs. 
These observations confirm that flares are detectable in 
the infrared portion of M dwarf spectra, which is much 
brighter in quiescence than the bluer portions of M dwarf 
spectra which are typically used to detect flares. 

Using a hotter chromosphere than previous one- 
dimensional static flare mo dels (e.g., iChristian et al.l 
I2003t iFuhrmeister et al.ll2010D . the ratios of Ca II K, He 



I A447lA, US, the Paschen lines, and Br 7 to H 7 can be 
relatively well reproduced. The generation of Ca II K 
in the hot, upper portion of the atmosphere is distinct 
from previous results, and is necessary to produce more 
Ca II K than Paschen series emission, which is observed 
during our strongest flare. This result confirms that in- 
frared emission is a useful constraint on the atmospheric 
heating during M dwarf atmospheres. 

The strength of emission from He I A10830A is not pre- 
dicted from our one-dimensional model, but including a 
detailed treatment of backwarming from the corona (e.g., 
lAllred et alll2006h may be warranted, based on solar re- 
sults. Modeling He I A10830A is also complicated by 
its different emission strengths compared to P/3 and P 7 
in the two flares on EV Lac, but these differences show 
that He I A10830A has potential to constrain different 
backwarming scenarios during a variety of flares. The 
time-evolution of the largest flare is not reproduced by 
our one-dimensional models, but a combinati on of mul- 
tiple models with different filling factors (e.g. IWalkowiczi 
120081 : iKowalski et al.ll2010T) or detailed radiative hydro- 
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Fig. 9. — The atmospheric temperature structure as a function 
of log(col mass) for the best fit model (log(col mass)y m i„ = 0, 
log(col mass)tr = —3.5, Ttr = 25 000K). The range of formation 
depth is shown for each of the emission lines; the lines are grouped 
by type and offset in Temperature for clarity. 



dynamic mqdelin 
lAllred et al.ll200~ ' 
flare emission. 



with non-thermal beam heating (e.g. 
may provide a better match to the 
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